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Surface-enhanced Raman spectroscopy (SERS) allows single-molecule detection due to 
the strong field localization occurring at sharp bends or kinks of the metal-vacuum inter- 
face. An important question concerns the limits of the signal enhancement that can be 
achieved via a judicious design of the surface. By using a specific example of a technologi- 
cally realizable nanopatterned surface, we demonstrate that while very high enhancement 
factors (~ 10 12 ) can be found for an ideal surface, these are unlikely to be achieved in 
laboratory samples, because even a minute, inevitable rounding-off strongly suppresses 
the enhancement, as well as shifts the optimal frequency. Our simulations indicate that 
the geometric enhancement factors are unlikely to exceed ~ 10 8 for real samples, and that 
it is necessary to consider the geometric uncertainty to reliably predict the frequency for 
maximum enhancement. 
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1 Introduction 

Raman spectroscopy is a powerful technique which has become one of the optical 
workhorses in analytical chemistry where it is used to detect and quantify a variety of 
chemical substances including biomolecules absorbed to a substrate. The 1974 report 
by Fleischmann et al. [1] on surface-enhanced Raman spectroscopy (SERS) created an 
avalanche in surface-enhanced spectroscopies [2]. Localized surface plasmons may dra- 
matically enhance the otherwise intrinsically weak Raman scattering [see panel (a) in 
Fig. [T] and this has paved the way for chemical detection in very dilute solutions, reach- 
ing the limit of single-molecule detection [3]. Initially, there was a strong emphasis on 
disordered and rough metal surfaces [4] ; however during the past decade the use of highly 
engineered surface plasmonic structures has been envisioned [5]. The latter development 
has to a large extent been fueled by the rapid development of a variety of nanofabrica- 
tion technologies, including high-resolution electron-beam lithography on the true nano 
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scale. As a consequence, surface plasmon subwavelength optics has in general developed 
tremendously [6], and in particular SERS has been observed for a large variety of nanos- 
tructures, such as nanorods, nanorings, triangular nanoprisms, nanocubes, nanocrescents, 
and nanorice [7]. 

Field localization is the key to understanding and further enhancing SERS. Structures 
with pronounced bends in their surface d'V will tend to localize the field also in the plane of 
the surface [8]. Singular points, such as sharp kinks and corners, may lead to nonanalytic 
or even mathematically divergent electric fields. Real structures are smooth, of course, 
at least on some microscopic length scale. All fabrication processes have their inherent 
upper bounds on the sharpness that can be fabricated reproducibly. As an example, 
state-of-the-art groove waveguides have sharp features down to a 10 nm scale [9]. In an 
ensemble of nanostructures there will always be a size dispersion or structure-to-structure 
fluctuations in some details of the geometry, for example the nanoscale rounding 5R [10] 
will differ from structure to structure, with possibly far-reaching effects on the SERS. 
Fluctuations between repeated features and variations in the smoothness of individual 
features have previously been addressed for fractal structures in a number of papers by 
Sanchez-Gil et ai, see e.g. Refs. [12,13]. The pronounced sample-to-sample fluctuations 
are potentially a serious challenge to the development of a technology for reliable use in 
quantitative sensing applications and bio-chemical analysis. 

In this paper we address another crucial issue: how are the high enhancement factors 
found for highly ideal surface structures affected by round-off effects? By means of nu- 
merical examples we shall illustrate the extreme precision needed in the fabrication of 
nanostructured surfaces, if one aims at SERS in the excess of ten orders of magnitude. 
Specifically, we consider the simple model geometry in panel (b) of Fig. [TJ Historically, 
this was the first structure considered theoretically to explain the large enhancement oc- 
curring for rough metal surfaces [4]. Today, we may easily envision realization of groove 
structures [14] closely resembling the geometry considered in panel (c), though state-of- 
the-art structures may still be somewhat larger in scale. As mentioned above, the points 
of intersection of neighboring cylinders may lead to a singularity in the electric field, and 
consequently arbitrarily large enhancement factors [see Eq. (0Q) below]. 



In the context of SERS it is customary to estimate the surface-enhancement factor 7 



where r is the position of the molecule. To a first approximation the surface-enhancement 
factor is independent on the intrinsic properties of the molecule. Here, d'V is the surface 
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bounding the volume "V occupied by metal, Ej is the incident field, and the total field is 
E t = Ej + E s with E s being the scattered component that arises in the presence of the 
metal surface. The surface-averaged enhancement factor applies to the case of a spatially 
uniform surface- coverage of molecules. 

Our theoretical study takes the wave equation for the electrical field as a starting point 

L) 2 

V x V x E(r) = — e(r, w)E(r) (2) 




where u is the angular frequency and c is the vacuum velocity of light. For structures as 
in Fig. [T] the spatially dependent relative dielectric function is of the form 



(3) 



where e(u) is the dielectric function of the metal. To make our simulations mimic a real 
material we use the bulk dielectric constant of silver (thus neglecting spatial dispersion 
near the metal surface), illustrated in Fig. [2] which shows standard handbook data for 
e(u>) [15] along with a numerical fit to the Drude expression 

e(u) = eoo + 4 

UJ Z + IUJ/T 

with ujp being the plasma frequency and r _1 the damping rate. This simple model accounts 
very well for the real part, however some deviations are seen in the visible for the imaginary 
part. The spectral position of resonances will thus be accurately predicted by using the 
Drude fit while we shall see that even small deviations in the imaginary part of e may 
influence the enhancement factor, Eq. (CQ), by orders of magnitude. 



3 Numerical results 

In the following we use a finite-element approach (Comsol Multiphysics) to solve Eq. (J2J) 
in a unit cell resembling panel (c) in Fig. [U We apply periodic boundary conditions in the 
direction parallel to the surface. Deep inside the metal and far out in the free space we 
employ perfectly-matched layers in the perpendicular direction [16], so that the domain 
has an extension 3> A, i.e. significantly exceeding the free-space wavelength. Eq. (j2J) is 
solved for a normal incident electrical field polarized perpendicular to the groove axis, see 
panel (b) in Fig. [H which strongly polarizes the surface by excitation of a surface-plasmon 
resonance. With the resulting scattered field, and hence E t , we subsequently evaluate 
Eq. (CQ) by a boundary integration routine. Our meshing is chosen to ensure convergence 
of the enhancement factor; this requires a particularly fine mesh in the groove region. As 
an example, resolving the peak in 7 at resonance has typically required a sub-nanometer 
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meshing. On a standard-equipped PC the calculation of a single near-resonance frequency 
point typically takes around 30 seconds. 

Figure [3] shows the wavelength dependence of the surface-averaged enhancement factor 
in the case of R = 15 nm and 5R = 0.1 nm. The solid line shows results based on 
the Drude model, Eq. (j3J), while the points are obtained from a simulation using the 
discrete frequencies where handbook data is available, see Fig. [2j First, we note that 
these geometrical parameters easily lead to surface-averaged enhancement factors in excess 
of 10 12 , i.e. twelve orders of magnitude, in the visible. Second, when inspecting the 
peak around A ~ 500 nm and comparing the two simulations, it is also clear that the 
Drude model, Eq. (jl]), overestimates the enhancement by perhaps up to three orders of 
magnitude. This overestimation may be traced back to the imaginary part of e shown 
in Fig. [2] where the Drude model seriously underestimates the broadening of the surface- 
plasmon resonance. Finally, when approaching the ultra-violet, resonances predicted by 
the Drude model are barely visible when employing the handbook data directly. While 
general trends can be extracted from a Drude description, these results emphasize the 
need for an accurate material-dispersion law for even semi-quantitative predictions of 
SERS. In particular, the surface-plasmon damping is crucial in the vicinity of resonances 
as clearly seen in the dashed traces where we for illustrative purposes have increased the 
damping rate r" 1 by a factor of two, five, and ten. 

While previous studies [4] have reported a spectrally broadened enhancement, our detailed 
finite-element simulations suggest the existence of highly peaked enhancement factors 
associated with the formation of a surface-plasmon resonance strongly localized in the 
groove. A related transverse localization has recently been reported experimentally for 
groove waveguides [14] and subsequently confirmed numerically [8]. Obviously, spatial 
localization is the route to high enhancement and especially one should minimize the 
overlap with the lossy metal causing Ohmic broadening of resonances. The nanoscale 
cutoff 5R is the geometrical knob that controls the strength of localization and the overlap 
with the metal. A stronger localization (5R — > 0) will usually be expected to squeeze up 
the resonance in frequency space. However, when increasing SR we actually observe a 
blue-shift which can be understood within the quasi-static framework of Ref. [17] [see 
their Eq. (4)] where an increase in metal volume V causes modes at uo to blue-shift to 
a frequency u> + Au with a less negative dielectric function, i.e. e(u + Au) > e(uj). In 
Fig. H] we illustrate this point for varying values of 5R at fixed R = 15 nm. For simplicity 
the results are based on the Drude model. Our results show that even for a relatively 
conservative value of the ratio 5R/R ~ 3% the enhancement has already dropped to 
7 ~ 10 8 compared to 7 ~ 10 12 in the case of 5R/R < 1%. 

Another geometric factor affecting the field enhancement is the radius of cylinders. In 
Fig. [5] we illustrate this point for varying values of R at fixed 5R = 1 nm, again using the 
Drude model. Our results show a red-shift of the resonance when increasing the radius 
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while the magnitude decreases slightly. The red-shift may be understood qualitatively 
from the scale invariance of Eq. ([2} in the absence of material dispersion. Figs. H] and [5] 
together illustrate how the enhancement is influenced by both length scales R and 5R in 
a nontrivial interplay. 

4 Conclusion 

In conclusion, we have presented finite-element simulations of the SERS enhancement 
factor for an example of a nanostructured surface, which represents the state-of-the-art 
technology. We find that a well-defined structure, with mathematically sharp bound- 
aries [see Eq. (J3J)], may exhibit huge enhancement factors, 7 ~ 10 12 . But any geometric 
smoothening rapidly decreases the enhancement, and even though a computational opti- 
mization may be possible by a painstaking survey of the parameter space (given that the 
dielectric function of the material is known very accurately - no small feat by itself!), it 
is doubtful whether the extremely high theoretical enhancement factors can be reached 
in practice and in a reproducible manner. 
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Figure 1: (color online) (a) Energy diagram for the Raman scattering process involving 
the absorption of an incident photon with energy ftw, and a subsequent emission of a 
photon with energy hw' via an intermediate virtual state in a molecule. The scattering is 
inelastic and the scattered light is Stokes shifted by Au = uj — uj' . (b) Metallic periodic 
structure composed of infinite half-cylinders of radius R. The normal incident electrical 
field Ej is polarized transverse to the cylinders, (c) Magnification of the groove structure 
indicating the nanoscale cutoff modelled by a finite radius of curvature 5R. 
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Figure 2: (color online) Dielectric function of silver. The data points show handbook data 
from Ref. [15] while the solid lines are fits to the Drude model in Eq. (j3J). 
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Figure 3: (color online) Wavelength dependence of the surface- averaged enhancement 
factor for a structure with R — 15 nm and 5R = 0.1 nm. The data points result from 
direct simulations at discrete frequencies using handbook data for silver, see Fig. El while 
the solid line shows the corresponding results based on the Drude model, Eq. (j4j). The 
dashed lines illustrate the suppression of the SERS effect for increasing damping rates 
r -1 . The inset shows 7 at the main peak where the electrical field is strongly localized in 
the groove. The electrical field points in a direction parallel to the superimposed electric 
field lines. 
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Figure 4: (color online) Wavelength dependence of the surface- averaged enhancement 
factor for a structure with R = 15 nm and varying nanoscale cutoff 5R. 
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Figure 5: (color online) Wavelength dependence of the surface- averaged enhancement 
factor for a structures with varying radius R of the semi cylinders and a nanoscale cutoff 
5R = 1 nm. 
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